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Abstract 

We study a version of the recently proposed modified F{R) Hofava-Lifshitz grav- 
ity that abandons the projectabiUty condition of the lapse variable. We discovered 
that the projectable version of this theory has a consistent Hamiltonian structure, 
T— I \ and that the theory has interesting cosmological solutions which can describe the 

' eras of accelerated expansion of the universe in a unified manner. The usual Hofava- 

Lifshitz gravity is a special case of our theory. Hamiltonian analysis of the non- 
projectable theory, however, shows that this theory has serious problems. These 
^0 ' problems are compared with those found in the original Hofava-Lifshitz gravity. A 

QQ . general observation on the structure of the Poisson bracket of Hamiltonian con- 

\ straints in all theories of the Hofava-Lifshitz type is made: in the resulting tertiary 

r~| ' constraint the highest order spatial derivative of the lapse N is always of uneven 

"V* ' order. Since the vanishing of the lapse [N = 0) is required by the preservation of the 

Hamiltonian constraints under time evolution, we conclude that the non-projectable 
r~| \ version of the theory is physically inconsistent. 

^ 1 Introduction 
in : 

^ ! Last year the so-called Horava-Lifshitz theory of gravity was proposed [1] (see also [21 E]). 

\ This theory is a candidate for a quantum field theory of gravity that aims to provide an 
^ ' ultraviolet (UV) completion of General Relativity (GR). At short distances it describes 
interacting nonrelativistic gravitons. Hofava-Lifshitz gravity exhibits anisotropic scaling 
of space and time coordinates 
> ! X ^ bx , t bH (1.1) 

^ ' with a dynamic critical exponent z = 1, 2, 3, . . .. In the UV regime the value of the critical 
■ exponent z is chosen so that the gravitational coupling constant is dimensionless. In 
[D + l)-dimensional space-time we have the scaling dimension [k^] = z — D. Thus the 
choice z = D is argued to ensure that the theory is power-counting renormalizable. For 
the usual case of 3-dimensional space, D = 3, we choose z = 3. 

The space-time manifold Ai is assumed to possess a foliation structure that enables 
one to define as a union of space-like hypersurfaces of constant time t. Due to 
the foliation the space-time is invariant under the foliation-preserving diffeomorphisms, 
whose infinitesimal generators are of the form 

Sx = at,x), St = fit), (1.2) 
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instead of the full diffeomorphism invariance of GR. For simplicity the topological struc- 
ture of space-time is assumed to be such that every leaf of the foliation is equivalent 
to a fixed manifold S: A4 = M x S. The preferred foliation of A4 enables the inclu- 
sion of spatial covariant derivatives into the action, which improve the UV behaviour, 
while avoiding time derivatives higher than the second order, which are known to produce 
problematic ghosts. 

At low energies and large distances the critical exponent is expected to flow to z = 1, 
so that the theory can coincide with GR. The Lorentz symmetry emerges at low energies 
as an accidental or approximate symmetry, but it is absent in the fundamental description. 

The theory comes in two flavors, with or without the projectability condition that 
requires the lapse to depend only on the time coordinate, = N(t). The projectability 
condition is one of the features that makes the theory differ from GR. Note, however, that 
many solutions of GR, as well as of non-projectable Hofava-Lifshitz gravity, respect the 
condition N = N{t) even thought the general theory does not. 

In the original theory an additional symmetry, the condition of detailed balance, is 
assumed. It defines the potential part of the gravitational action in terms of a variation 
of a D-dimensional action on the spatial hypersurface with respect to the spatial metric. 
The purpose of the detailed balance condition is to reduce the number of independent 
couplings to 3 — otherwise there are 9 independent couphngs. We, however, do not 
assume this condition, so all terms that have appropriate scaling properties and that are 
covariant under foliation-preserving diffeomorphisms can be included. 

This theory has received a lot of attention and many potential problems have been 
discovered. Some of the problems are very serious. First it was found that GR is not 
recovered at large distances if the detailed balance condition is assumed [H [5]. A "phe- 
nomenologically viable" version of the theory without the detailed balance condition was 
soon introduced [6]. Due to the reduced diffeomorphism symmetry group there is an ad- 
ditional "half" scalar degree of freedom in Hofava-Lifshitz gravity. It has been shown to 
be strongly coupled at all scales by considering perturbations about a reasonable vacuum 
[7], regardless whether the detailed balance is assumed or not. This suggests that per- 
turbative GR cannot be reproduced in Hofava-Lifshitz gravity [7j, and that the theory 
could be ruled out by existing observations on the gravitational radiation of binary pul- 
sars, which agree with linearized GR. The low-energy regime of the theory was further 
analyzed in Ref. [H [9] where problems with instability and strong coupling of the extra 
degree of freedom were found. Since then these problems have been confirmed in various 
papers. Although the non-projectable version of the theory may not give GR as the limit 
at large distances, some other scenarios, such as the chameleon, may solve this problem. 
A "healthy extension" of Hofava-Lifshitz gravity was proposed in Ref. [10] that is argued 
to be free from at least some of the pathologies of the original theory, since the extra 
scalar mode has a healthy quadratic action. This is achieved by adding terms that involve 
the spatial 3-vector N~^'ViN into the action. We assume that the Hamiltonian takes 
the canonical form [see (12.191) ] which rules out such terms. Hamiltonian formalism of the 
healthy extension has been studied in Ref. [TIJ . 

Renormalizability of the Hofava-Lifshitz gravity has been investigated beyond the 
power-counting scheme in Ref. [12j. 

When the projectability condition is assumed, the theory has a quite simple and con- 
sistent Hamiltonian structure. The algebra of constraints was shown to be closed for 
2; = 1,2 in Ref. [2], and this holds for a higher scaling exponent z as well. The Hamil- 
tonian structure of Hofava-Lifshitz gravity without the projectability condition has been 
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analyzed particularly in Ref. [131 [H]. The non-projectable theory is physically inconsis- 
tent for generic couplings [H], including the case with detailed balance p!3]- In the case 
of low-energy effective action a consistent set of constraints can be obtained by imposing 
an additional constraint (tt = 0) [HI [151 [16] - 

Recently we proposed the modified F{R) Horava-Lifshitz gravity [TTlllH] that combines 
the interesting cosmological aspects of F{R) gravity and the possible UV finiteness of 
Hofava-Lifshitz gravity. In particular, we demonstrated that the solution of spatially-flat 
FRW equation has two branches: one that coincides with the usual F{R) gravity for a 
certain choice of parameters, and one that is totally new and typical only for Horava- 
Lifshitz gravity. It was shown that unlike to standard Hofava-Lifshitz gravity, our F{R) 
Hofava-Lifshitz gravity enables the possibility to unify the early-time inflation with the 
late-time acceleration in accord with the scenario of Ref. [19]. In this paper we present 
the Hamiltonian analysis of the non-projectable version of this theory, where the lapse 

depends also on the spatial coordinates: N = N{t,x). Expectedly the Hamiltonian 
structure of this theory turns out to be more complex than in the projectable case. Our 
analysis should be of interest to everyone interested in the Hamiltonian formalism of 
gravity, and of modified gravity in particular. 



2 Hamiltonian analysis 
2.1 Action 

We assume the ADM decomposition of space-time [20] (for reviews and mathematical 
background, see [21]). The metric tensor of space-time is decomposed in terms of the 
ADM variables as 

V^dx'^dx^ = -{N^ - NiN')dt^ + N,{dtdx' + dx'dt) + gijdx'dx^ , (2.1) 

where N is the lapse, N'^ is the shift vector, gij is the spatial metric tensor, and x\ i = 1, 2, 3 
are spatial coordinates on the t = constant hypersurface S^. The covariant derivatives 
defined by the metric tensors ^(7^,^ and g^j are denoted by V^'^'* and Vj, respectively. The 
extrinsic curvature of the hypersurface is 

K^J = ^ {9^J - 2V(,,Ar,)) , (2.2) 

where the dot denotes the derivative with respect to time t. The scalar associated to the 
extrinsic curvature is denoted hj K = g^^Kij. The (intrinsic) curvature of the space Sf 
is defined by the spatial metric gij in the usual manner. The natural invariant volume 
element of space-time is decomposed 

d^x = dtd^x^N . (2.3) 

The action of the non-projectable version of the modified F{R) Hofava-Lifshitz gravity 
is defined similarly as in the projectable case [T7] : 

SF = ^j dtd'x^NF (^rJ , 
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Here A and are constants, n'^ is the unit normal to the spatial hypersurfaces Sj, and 
Cpiigij) is a function of the three-dimensional metric gij and the covariant derivatives Vj 
defined by this metric. The crucial difference compared to the theory we proposed and 
analyzed in Ref. [17] is that the lapse function does not obey the projectability condi- 
tion, i.e. it depends also on the spatial coordinates, = N{t,x). For the Hamiltonian 
analysis of modified Hofava-Lifshitz-like F{R) gravity, which is a special case of the gen- 
eral projectable theory ^7\, one can see Ref. [HI 122]. This special case with the further 
restriction to the parameter value /i = has been proposed and analyzed in Ref. [23] . 
Yet another special case, given by F{^R) = *R and Cn^gij) = —f{R), has been studied in 

Ref. m- 

By introducing two auxiliary fields A and B we can write the action (12.41) into a form 
that is linear in '^R: 

Sp = ^ j dtd^x^N B (^R - ^) + F{A) . (2.5) 

The variation with respect to B yields ^R = A, which can be inserted back into the action 
(12. 5p in order to produce the original action (12. 4p . The variation with respect to A yields 
B = F'{A), where F' denotes the derivative of F with respect to its argument. Thus (12. 5 p 



reduces to the action fl2.4p when these equations of motion are imposed. 

First we rewrite ^i? in (12. 5p in a more explicit and useful form (see (12. 4p for the definition 
of ^R). The unit normal to the hypersurface in space-time can be written in terms of 

the lapse and the shift vector as = (n°,?T,*) = (^jf,~^^- The corresponding one- form 

is 71^ = —NV^i^H = {—N, 0, 0, 0). The term in (12. 4p that involves the unit normal can be 
written as 

V^f) {n^V^^'^n'' - n'^V'^^^n^) = V^f^ (n^K) - j^AN , (2.6) 
where the spatial Laplacian is A = g^^ViVj. Thus we can rewrite '^R 

^R = K,,g^^''K,i + 2/iVf (n^K) - ^AN - Cnig.j) , (2.7) 
where the "generalized De Witt metric" was introduced 

Introducing (12.70 into (12.51) and performing integrations by parts yields the action 

Sf = -^j dtd^x^^N [B {Ki,g'^''Kki - CR{g,,) - A) + F{A)\ 

-2nK (b - N'diB^ - 2i^NAB^ , (2.9) 

where the integral is taken over the union U of the t = constant hypersurfaces with t 
the boundary integrals over dU and d'Et vanish. 



over some interval in R, and we have written Nn'^Wfj, B = B ~ N^diB. We assume that 



2.2 Hamiltonian and momentum constraints 

In the Hamiltonian formalism the field variables g^j, N, Ni, A and B have the canonically 
conjugated momenta vr*-', ttjv, tt', tta and tib, respectively. Because the action does not 
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depend on the time derivative of A^, A^* or A, their conjugated momenta are the primary 
constraints: 

7r^(a;)^0, 7r^(a;)^0, 7r^(a;)^0. (2.10) 
For the spatial metric and the field B we have the momenta 

6Sf 



7b ~ 



[bG'^'^Ku - ^g'' [is - N'd,B^] , (2.11) 



VTB = -^ = -^^K. (2.12) 



We assume 7^ so that the momentum (12.121) does not vanish. When A 7^ 1/3, the 
generalized De Witt metric (12. 8p has the inverse 

Gijki = - {9ik9ji + 9ii9jk) - ^ 9ij9ki , GijkiQ''^"''' = ^'(T^'j) ■ (2-13) 

However, as long as /i 7^ 0, the invertibility of (12. 8 p is not that significant in our theory, 
because K is given by (12.121) . K = —^^j-jf^B-, and therefore we have 

Q'^^'K^, = K^^ + ^^9'''nB . (2.14) 

The case /i = will be discussed later in Sec. 12.51 

The Poisson brackets are postulated in the form (equal time t is understood) 

{g,,{x),7,>^\y)} = 5^;;5l5{x-y), 

{N{x), TXN{y)} = 5{x - y) , {N,{x), Ti^{y)} = 5l5{x - y) , 

{A{x), ^A{y)} = 6{x - y) , {B{x), TXB{y)} = 5{x - y) . (2.15) 

All the other Poisson brackets between the variables vanish. We shall continue to omit 
the argument {x) of the fields when there is no risk of confusion. 

In the following analysis we may lower and raise spatial indices with the 

spatial metric gij and its inverse e.g. Hij = gikgjiir^'', and we will denote 

7T = g^j7T''. (2.16) 

For some integrals over the space S^, we will routinely perform integration by parts, 
then apply the divergence form of the Stokes theorem and assume that the resulting 
boundary integrals over dJ^t vanish so that they can be ignored. This can be justified 
by assuming appropriate boundary conditions for the variables (asymptotic behaviour at 
infinity), similarly as in GR. 

In order to obtain the Hamiltonian, we first solve (I2.11l) - (l2.12p for Kij and B, 

K., = — (^vr., - -,.,.rj - -g.,.s\ , (2.17) 
B ^ N^^^B-N^^{^^^.^l^B.:). (2.18) 

and further obtain cjij = 2NKij + 2V(^iNjy Therefore both gij and B are dynamical 
variables and no more primary constraints are needed. The Hamiltonian is then defined 

H = J d^x (n'^gij + tcbB^ -L = j d?x {NUq + NiW) , (2.19) 



where the Lagrangian L is given by the action (12. 9p . Sp = J dtL, and the so-called 
Hamiltonian constraint and the momentum constraints are found to be 
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+ ^ [B {CR{g,,) + A)- F{A) + 2/iAi?] , 

= -2d,7r'^ -g'^2dk9ji-d,gH)7r''' + g''d,B7rB, (2.20) 
respectively. We define the total Hamiltonian by 

Ht = H + j d^x {Xnttn + AiTT* + Xatta) , (2.21) 

where the primary constraints (I2.10p are multiplied by the Lagrange multipliers Xn, Aj, 
A^. The total Hamiltonian (I2.2ip generates the time evolution of dynamical variables: 

f{x) = {fix),HT}. (2.22) 

The primary constraints (I2.10p have to be preserved under the time evolution of the 
system: 

ttn = {vtat, Ht} = -"Ho , 
r = {tt\Ht} = -W, 

TTA = {'Ka,Ht] = ^N{-B + F'{A)) . (2.23) 

Therefore we impose the secondary constraints: 

'Ho{x)^Q, W{x)^0, 

^a{x) = B{x) - F'{A{x)) ^0. (2.24) 

The Hamiltonian constraint Tioix), the momentum constraints W{x) and the constraint 
^a{x), are all local. It is convenient to introduce globalized versions of the Hamiltonian 
and momentum constraints: 



$o(r?) = J d'xr]no^O, 

$5(ei) = [ d'xm'^0, (2.25) 



where rj and ^j, z = 1, 2, 3 are arbitrary smearing functions that vanish rapidly enough at 
infinity — the choices t] = S{x — y) and = SjS{x — y) will produce the local constraints 
Hq and 1-L\ which in turn imply the smeared constraints. 

2.3 Consistency of the secondary constraints under dynamics 

The total Hamiltonian fl2.2ip can be rewritten in terms of the Hamiltonian and momentum 
constraints (12.250 as 

Ht = %{N) + + j d^x (Aatttjv + X.y + A^tt^) . (2.26) 



The consistency of the system requires that also the secondary constraints $o(^), ^sid) 
and ^a{^) have to be preserved under time evolution generated by the total Hamiltonian 

$5(6) = {$5te),$o(iV)} + {$5te),<^>5(iVi)}~0, (2.27) 

where we have used the fact that the constraints ttjv and vr* have strongly vanishing Poisson 
brackets with every constraint, and that 

{$5te),vrA} = 0. (2.28) 

We need to calculate the rest of the algebra of the constraints under the Poisson bracket. 
The Poisson brackets between the momentum constraint ^s{^i) and the canonical vari- 
ables are 

{$5(6), 5} = -CdiB, 

{^s{ik)y'} = -dk{e^'') + 7r^''dke + 7i"'dke, (2.29) 

and trivially zero for A and tta, 

{$5(6), ^} = , {$ 5(e.), tta} = . (2.30) 

Thus ^si^i) generates the spatial diffeomorphisms for the variables 5, tt^, (jfjj, tt*-', and 
consequently for any function or functional constructed from these variables, and treats 
the variables A, t^a as constants. By using this result fl2.29p - fl2.30l) we obtain the Lie 
algebra of the generators ^si^i)'- 

{$5(6), <^siv^)} = $s(e5,^i - V'dj^^) ^ . (2.31) 
Then we calculate their Poisson brackets with the Hamiltonian constraint $o('?)- 

{<!>sm, %{V)} = %{ediV) ~ , (2.32) 

which tells us that Hq is a scalar density under spatial diffeomorphism. The momentum 
constraints W{x) will be first-class everywhere and consistent under time evolution. 
Then consider the constraint ^a{x), whose nonvanishing Poisson brackets are 

{$.(.), $o(^)} = -^;^(^- + ^^-- 

{<^Aix),<l>sm = CdiB, 

{^A{x),My)} = -F"{A{x))5{x-y). (2.33) 

Thus, in order to satisfy the consistency conditions fl2.27p . we have to impose the tertiary 
constraint 

N^d.B - (^TT + ^—^Bttb] - XaF"{A) ^ . (2.34) 
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Since F"{A) = would essentially reproduce the usual non-projectable Hofava-Lifsliitz 
gravity, we assume that F"{A) ^ 0. The first two terms in (12.341) . i.e. the expression for B 
in (I2.18p . does not vanish due to the established constraints (I2.10p and (I2.24p . Therefore 
(I2.34P is a restriction on the Lagrange multiplier Aa, and we can solve it: 



A, 



( N%B - N— (—7^ + - 
F"{A) \ ^\3fi 



3A 



6/i2 



-Btvb 



(2.35) 



Introducing fl2.35p into the Hamiltonian f l2.26p ensures that the constraint ^a{x) is con- 
sistent. 

Finally we consider the Hamiltonian constraint ^o{v)- The Poisson bracket with the 
primary constraint tta{x) vanishes in (I2.27P 



(2.36) 



Also the Poisson bracket of the Hamiltonian constraint $o(^) with itself has to be calcu- 
lated in order to check its consistency. Since the Poisson bracket is antisymmetric with 
respect to its arguments and the Hamiltonian constraint does not depend on the primary 
constraints (I2.10p . we can write the Poisson bracket of Hamiltonian constraints as 



or 



{$o(0,*o(r/)}= / d'x 



+ 



Sgij{x) 6'n-''-j{x) 6B{x) dnsix 



(e ^ v) (2.37) 



+ {$0(0, nBix)}{Bix), %ir])}) - ^ v) ■ (2.38) 



Further simplification follows from the fact that the Poisson brackets of the Hamiltonian 
constraint ^o{v) with the field variables Qij and B are proportional to rj, i.e. of the form 
rif where / is a function of the canonical variables: 

^2 



{B{x),%{r^)} 



( 1 

—ri — TT 

^f9 V3/^ 



1 -3A 



3/i 



2r]K, 



BllB 



(2.39) 
(2.40) 



Therefore, the parts of the Poisson brackets of the Hamiltonian constraint $0(0 with 
the momenta vr*-' and tt^ that are proportional to ^ will contribute to terms that are 
proportional to ^r] = t]C,. These terms will necessarily cancel out due to the antisymmetry 
of fl2.38p under ^ O 77. Thus, in order to calculate f l2.38p . we only need those parts of the 
Poisson brackets with the momenta that contain spatial derivatives of ^. These Poisson 
brackets are found to be 



{$o(0,^*'(^)} 



K 



9' 



V9 



1 

2B 



TTfc/VT - -TT 



1 

6/i 



1 -3A 
24/i2 



Bn^ 



B 



1 



-B{CR{gu) + A)--F{A)+ii/\B 



V9 VB 



jk 



1 

3' 



1 

3/i 



— tc^^tcb 



(2.41) 



+ 1 y d'y^^y) [B{y) 



6Cnigkiiy)) ^ n^ ^^Bjy] 
6gij{x) 6g. 



X 
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and 



V9 



1 

52 



TCijTT 



12/i 



3A 2 



where the last term was obtained from 
Ik k 



(2.42) 



(2.43) 



We find immediately that only the last term in both (12.411) and fl2.42p give nonzero 
contributions to the Poisson bracket fl2.38p . 

Then we must calculate the variations of the potential part Cfjlgij) and AB with 
respect to the spatial metric in (I2.4ip . For this we need the variations of the involved 
geometrtic quanties. The variations of the connection coefficients T^j, the Ricci tensor Rij 
and the scalar curvature R are given by 



6R 



6R 



1 1 
-g'"^ {VkVi^gij + VkVj^gii - VNj^gki) - ^^5gij, 



-R^g,, + W^g,,-g''ASg,,. 
For AB we obtain the variation 



(2.44) 



6AB 



,W'B + -g'^ {y.,5gu, + V j5giu - Vu^gi,) V'B 



and then in fl2.4ip we have 



5AB{y) 

Sgijix) 



2fi^ [-2^V'^'V^^B - V^'^V^^B 



+ ^g''{Vk^V'B + ^AB) 



(2.45) 



(2.46) 



where aij represents any tensor, and the first term and the last term will be cancelled 
for the familiar reason. Thus we obtain the general result for the Poisson bracket of 
Hamiltonian constraints fl2.38p : 



2 
B 



3yU 



gijT^B 



\.^'-^b.:)a,]-^,^,). 



6C 



nygki 



[y)) 



6gij{x) 



(2.47) 



Further progress can be obtained by specifying the form of the potential part C[i{gij) in 
the action. 

For Cji{gij) let us start with the simplest case that could be the low-energy effective 
potential 

Cnig.j) = ao + aiR , (2.48) 
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where ao and ai are coupling constants. f l2.48p gives in (12 as a "deformed scalar 
curvature + constant": 



= K'^Kij - - aiR + 2^iVf (n'^K) - '^AN - ao 



(2.49) 



When A = l,yU=l,ai = — 1 and ao = 0, this '^R reduces to the scalar curvature '^R of 
space-time. For asymptotically flat spaces we assume ao = 0. For (I2.48P we obtain 



Sgij{x) 



(2.50) 

where Ojj represents any tensor, and again the parts that contain no derivatives of ^ will 
cancel out in the Poisson bracket (I2.47p . This expression contains first and second order 
spatial derivatives of ^. Thus we obtain the result 



where we have defined: 



(2.51) 



CI 



2ai7r 



9 



2(ai + l) l-2ai-3A 



-vr 



3/i 



Btvb 



4(ai - /i)7r*^' + g'^ 



o ^ H o B'^B 

6 6H 



B 



(2.52) 



The tensor C2 is symmetric and so is the part of C| that multiplies Vj-B in the definition 
(I2.52p . Note that in the case ai = — 1,A = 1 the coefficient of the A^-term vanishes in 
(I2.5ip so that C2 — 2ai7r*-'. Finally we can write (I2.5ip into a form where the integrand 
is proportional to ^: 



{$o(0,$o(r7)} = jd^x^ (EiViV + Eov) 



(2.53) 



where we have defined 



El 

Eq 



2 {V,Cl^ - Ci) , 



(2.54) 



Note that the terms involving the second spatial derivative V(iVj)?7 cancel each other 
independently of the form of C2 . Now the condition that is necessary for the preservation 
of the local Hamiltonian constraint 1-Lq{x) can be obtained by inserting 1^ = 5{x — y) and 
7] = N into the Poisson bracket (12.53p . Thus, in order to ensure the consistency of the 
Hamiltonian constraint 'Ho{x), we must impose the constraint 

E{ViN + EoN ^ , (2.55) 

where we have factored out some parts of the E^s that vanish due to the momentum 
constraints by defining: 

Ei = E{ + 2aiW 



_ 2{a, + l) ^,^ + 4(-ai + /i) - 



B 



1 - 2ai - 3A 
3/.t 



Eq 



BV'tcb + 
Eo + aiViW = ^V,E[ . 



3A + /i — 2ai — 1 
3/1 



+ ai V'Bttb 



(2.56) 
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This is a homogeneous first-order partial differential equation for the lapse N. As such it 
always has the solution = 0. Due to the relation (12.541) of Eq and El, we can rewrite 
the constraint f l2.55p as a divergence (after multiplying it by 2A^) 

Vi (^N^Ei^ ^ . (2.57) 

Let us next consider the case with the critical exponent z = 3 that could provide an 
UV complete theory. For the potential Cnlgij) there are many terms that have the same 
scaling dimension as the kinetic terms under the anisotropic scaling (II. ip with z = 3. 
Such terms are, for example, the terms quadratic in curvature 

V'RViR, V'R^'^ViRjk, V'R^'^VjRki (2.58) 

and 

RAR, R'^ARij, R'^ViVjR, (2.59) 

which modify the propagator in addition to providing interactions, and the terms cubic 
in curvature that are pure interactions 

i^^ RR'^Rij, R)R\R\. (2.60) 

Many of the terms are related to each other due to the properties of the Riemann tensor 
in three dimensions (Weyl tensor vanishes), the Bianchi identity, and integration by parts. 
Indeed only two terms of the types (12.580 and (I2.59P are independent, which we choose 
to be 

V'R^^ViRjk, R^R- (2.61) 

The most general potential that contains all the independent renormalizable and super- 
renormalizable terms, while maintaining the canonical form of the Hamiltonian (I2.19p . 
is 

^R{9ij) = "0 + tti-R + "2-^^ + a^R'^Rij + a^R^ + a^RR'^Rij 

+ aQR'jR\R\ + a7V'R'''ViR,k + agRAR . (2.62) 

Under the variation of the spatial metric the terms involving second order spatial deriva- 
tives of curvature (I2.59P contain spatial derivatives of the variation 6gij up to fourth order. 
Likewise, the variations of the terms (I2.58P contain spatial derivatives of the variation Sgij 
up to third order, and the variation of (I2.60p up to second order. This means that the 
variation of the potential (12.620 contains spatial derivatives of the variation 6gij up to 
fourth order. Thus the Poisson bracket of Hamiltonian constraints has the form 

{$o(0, Mv)} = J d'xT] (cf V,,He + cf'^i.k^ + c'^v^,^ + c{v,i) -{i^v), 

(2.63) 

where C„'s are symmetric tensors consisting of the canonical variables and their spatial 
derivatives, and we denote Vij = V(iVj) etc. For arbitrary couplings am the tensors C„ 
have quite complicated forms, which we do not present here. After integration by parts 
we obtain 

{$o(0, Mv)} = f d'x^ [EfV^.k^l + E^^^JV + E\V,7^ + Eq//) , (2.64) 
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where we have defined 

Eq = VijkiCi^^ — ^ijkC-J^ + VjjCj"' — VjC*! . (2.65) 

The resulting tertiary constraint is of the form 

Ei^^VijkN + Ei^V.jN + E\ViN + E^N ^ . (2.66) 

The condition f l2.66p is again homogeneous in and is therefore satisfied by = 0. Note 
that the fourth order spatial derivative of cancels out similarly as the second order 
derivative did in the case fl2.48p . More generally for every even n, C„ contributes to Ej. 
for k < n, but never to En- For every uneven n, the contribution of Cn to En is — 2C„. 
Thus the highest spatial derivative of N in the constraints like f l2.55p and f l2.66p is always 
of uneven order. This is clearly a very general result that holds in any theory where 
the Poisson bracket of Hamiltonian constraints is of the form f l2.63p with nonvanishing 
coefficients Cn usually up to n = 2; + 1 — including the usual Hof ava-Lifshitz gravity. The 
constraints fl2.57p and f l2.66p will be discussed further in Sec. 12.71 

2.4 Elimination of the auxiliary field A 

According to the Poisson brackets between the constraints fl2.3ip - fl2.33p and fl2.36p . we 
can set the second-class constraints 7r^(a;) and ^a{x) to vanish strongly. Then the Hamil- 
tonian constraint $o(^) "will have weakly vanishing Poisson brackets with every constraint 
except with itself fl2.64p . The momentum constraints $s(^j) "will be first-class assuming 
its Poisson bracket with possible additional constraints are weakly vanishing too. To this 
end, we replace the Poisson bracket with the Dirac bracket, which is given by 

{f{x),h{y)}^^ = {f{x),hiy)} + / ^''^ F"{A{z)) i{fi^)^M^)}{'^A{z),h{y)} 

-{f{x),^A{z)}{7rA{z),h{y)}) , (2.67) 

where / and h are any functions of the canonical variables. Assuming we can solve the 
constraint $a(^k) = 0, i.e. B = E'{A), for A = A{B), where A is the inverse of the 
function E', we can eliminate the variables A and vr^. Thus the final variables of the 
system are gij, vr*-', B, hb- The lapse N and the shift vector Ni, together with Xn and Aj, 
are non-dynamical multipliers. Since every dynamical variable has a vanishing Poisson 
bracket with the constraint tta, the Dirac bracket fl2.67p reduces to the Poisson bracket, 

{f{x),h{y)}j,s = {f{x),h{y)}. (2.68) 

2.5 The cases F" = and ^ = 

For F" = 0, i.e. F{'^R) = ci'^R + cq with constants ci and cq, the action ([23]) with fl2:62|) 
reduces to the action of the usual Hofava-Lifshitz gravity without the detailed balance 
condition, when the couplings are redefined: Ci/k^ — >• I/k^, (ciOo — Co)/k^ — )■ and 
Cittn/ n"^ — )■ a^, for n = 1, 2, . . . , 8. The value of the parameter fi is irrelevant, since the term 
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involving (12. 6p is a total divergence and hence it can be dropped. For Hamiltonian analysis 
of non-projectable Hofava-Lifshitz gravity see Ref. [2] and references therein, recalling the 
observations about the general structure of the Poisson bracket of Hamiltonian constraints 
made in Sec. 12.31 

Let us then consider the case fi = 0, and assume F" 0. As we mentioned this is a 
generalization of the (second) model proposed in Ref. |23] without the detailed balance 
condition. The additional primary constraint 7Cb{x) ~ has to be introduced in addition 
to (I2.10p . The canonical momenta are 



1 

^" 2 



^BG^^'^^Kki . (2.69) 



The field B is non-dynamical. For A ^ 1/3, we obtain the Hamiltonian and the momentum 
constraints 



W = -2Vj7r^^ (2.70) 



•Ho = ^^7r'^g,jki7i'' + ^[BiCRig,,) + A)-FiA)] 



The preservation of the constraint 'Kb{x) under time evolution implies the secondary 
constraint 

The constraints Xai.'^) = {t^a^^a^t^b^^b) are second-class and the corresponding Dirac 
bracket is defined 



= -^^^r''G^Jkl7l'' - ^ {CR{g,,) + A) ^ . (2.71) 



{/(x), /^(?/)}db = {/N, Ky)} - j A'zA'z'ifix), Xa{z)}C'^\z, z'){x,{z'), h{y)} 

(2.72) 

where C"^(z, z') = C"'''{z)5{z — z') is the inverse of 

Cab{x, y) = {xaix), Xb{y)} = Cab{x)6{x - y) . 
The nonvanishing components of these antisymmetric matrices are 

Cu{x) = F"iA), C,,{x) = ^, 

C^six) = 1 , C,,{x) = -^n^^g.-^m'^i , (2.73) 



and 



C^\z) 
C'iz) 



2TT^igijkiTr''^F"{A) + gB^/K^ 



27i'^gijMrr'''F"{A)+gByK^ 

gB^K^ 
' 2n'i gi^kiTi'^'F" (A) +g By 



C^Uz) = ^9F"{A)BVn- 

^ ' 271^1 gi^kiT^'^'F" {A) +g By K^' ^ ' 
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Since the constraints xi = and Xs = have vanishing Poisson brackets with the 
canonical variables Qij and vr*-', and we have C'^'^{x) = 0, we again find that the Dirac 
bracket reduces to the Poisson bracket (12.681) for any functions of Qij and ir'^K The auxiliary 
fields can be eliminated from the action by solving ^a{^) = and ^si^) = for 

A = A (^^n'^g^^.in''^ ^CR{g,,)^ , B = F'(i) , (2.75) 

where A satisfies 

-^jr^,n-G,Mrr'' - ^ (Cnig.,) + a)=0. (2.76) 
Then the Hamiltonian constraint can be written in terms of Qij and vr*-': 

-Ho = ^^\~y 'Q^,kiri''' - ^F{A) . (2.77) 



^F'{A) 



The seemingly simple form of this Hamiltonian can be deceiving, since the form of the 
function A can be complex. For example, when F{R) is a polynomial of degree n, then 
A is a root of a polynomial p{A) of degree 2(n — 1) + 1. 

The momentum constraints fl2.70p clearly generate the spatial diffeomorphisms of the 
canonical variables, since "H* are identical to the momentum constraints of GR. It is also 
clear that the Poisson bracket of Hamiltonian constraints has the general form (12.631) when 
the potential takes the general canonical form (I2.62p . 

Since Q^^^^Qki = (l~3A)(y'*-', the value A = 1/3 implies the primary constraint 7r(a;) 
similarly as in the usual Hofava-Lifshitz gravity. Then the Hamiltonian and the momen- 
tum constraints are found to be as in (I2.70p but with i^'^^QijkiT^^^ replaced by iTijir'^K The 
preservation of ti{x) ~ under time evolution is a nontrivial matter that is not discussed 
here. 



2.6 On the class of theories with a projectable lapse 

For comparison we briefiy consider the class of theories with a projectable lapse N = N{t) 
(see [TTj for more details). In this case there is only one (global) Hamiltonian constraint 
$0 = / d^xTio = $o(l)- Now the Poisson bracket of Hamiltonian constraints vanishes, 

{<I>o,<l>o} = 0, (2.78) 

and there is no need for additional constraints like the one (I2.66P in the class of theories 
with a non-projectable lapse. 

As explained in Sec. 12. 4[ the second-class constraints ttaIx) and ^a{x) can be set to 
vanish strongly by utilizing the Dirac bracket (I2.67p . which reduces to the Poisson bracket 
(I2.68p . As a result the Hamiltonian constraint $0 and the momentum constraints ^s{^i) 
are first-class, since they have vanishing Dirac brackets with every constraint. Finally the 
total Hamiltonian is the sum of the first-class constraints 

Ht = N% + ^siNi) + Ajvvrjv + j d^ajA.vr' . (2.79) 

In the case /i = discussed in Sec. 12. 5[ the projectability condition has the exact same 
effect as in the case 7^ considered above. 
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2.7 Interpretation and analysis of the tertiary constraints — 
conditions on the lapse or not? 

Let us discuss the conditions f l2.57p and fl2.66p . In order to obtain real dynamics, the 
requirement 7^ is physically necessary, and if such a solution does not exist, the 
system would have to be concluded as physically inconsistent. Such inconsistency has 
been shown to exist in the usual non-projectable Horava-Lifshitz gravity with a generic 
potential |T1], where only the solution = exists when asymptotically flat spaces with 
appropriate boundary conditions are considered (see [U [131 El E| for other analyses) . When 
the lapse is fixed to = 0, the Hamiltonian constraints of these theories are generically 
second-class. The special cases of GR and ultralocal theory are known to be consistent 
and possess first-class Hamiltonian constraints. Other cases of the potential can at most 
have first-class Hamiltonian constraints in some parts of space, but not everywhere [26j . 
Imposing additional constraints may, at least in some cases, provide a way to define a 
consistent theory. For example, the case of low-energy effective potential of the form R 
provides a consistent theory when the additional constraint vr = is imposed [HI [151 HB] , 
and even to be equivalent to GR in the gauge vr = 0. Problems related to the modification 
of the Hamiltonian constraint of GR have also been studied in Ref . p7l |28] . 

For asymptotically fiat spaces we assume the following asymptotic behaviour of the 
canonical variables in asymptotically flat coordinates {C = constant): 



which are an extension of the standard boundary conditions of GR [25] that have been 
used in the usual Hof ava-Lifshitz gravity [H] . The auxiliary field A is expected to behave 
as ^i? for obvious reasons. 

The interpretation of the conditions fl2.57p and more generally fl2.66p is similar com- 
pared to the analysis of the usual Horava-Lifshitz gravity given in Ref. [T4j. They are 
generically conditions on the lapse N and the only solution that satisfies the appropriate 
boundary conditions in asymptotically flat space is = 0. Mathematically the constraint 
A^ ~ is perfectly acceptable, and setting the Lagrange multiplier = ensures that 
this constraint is preserved under dynamics. But as we have stated, this constraint ruins 
the chance of having any meaningful gravitational dynamics. A^ = implies that the 
Hamiltonian vanishes in the gauge A^j = 0, which means that every function is a con- 
stant of motion. Moreover, it is unlikely that the constraint f l2.66p could be satisfied by 
introducing additional constraints that are both consistent under dynamics and do not 
constrain the canonical variables too much. Such tertiary constraints, together with the 
existing primary and secondary constraints, would essentially have to imply En ^ in 
fl2.64p -( l2.66p . We will demonstrate this by considering the simple special case fl2.55p . 

In the low-energy effective case (I2.48P with only one nonvanishing coupling ai the 
condition fl2.57p cannot be satisfied by introducing a simple additional constraint like 
vr = in the usual theory [HI [15l [26l [16]. Instead we could try to introduce more 
complex constraints to serve the same purpose. In order to satisfy f l2.57p the following 




(2.80) 
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three constraints could be imposed: 




2(^1 + 1) 
3 

1 - 2ai - 3A 



3/i 




V'Btib^O. (2.81) 



It would require constraints like tt{x) ^ 0, 7rB{x) ~ and B{x) ^ constant (i.e. "pro- 
jectable B") to make the El vanish due to relatively simple constraints. But we can see 
from fl2.39p - fl2.42l) that the preservation of 7t{x) ^ and nBix) ~ under dynamics would 
impose additional (quartic) constraints on the canonical variables. These constraints most 
likely constrain the system too much, leaving too little physical degrees of freedom for 
gravitational dynamics. If we assume that only two quartic constraints would be needed 
and that all these five constraints would be second-class and no more constraints are 
needed to ensure their consistency, and that the variables A and tta can be eliminated 
(see Sec. 12. 4p . we find that the (maximum) number of physical degrees of freedom would 
be 3/2 at each point. This is one degree of freedom less than in the usual Hofava-Lifshitz 
gravity and one half degree of freedom less than in GR. Thus no more than three con- 
straints, e.g. fl2.8ip . can be introduced to satisfy the constraint (12.551) or fl2.57p . But the 

consistency of the three extra constraints (12.811) under dynamics is problematic — El is 
a long and complex expression that appears to imply further (quartic) constraints, which 
would again lead to the lack of degrees of freedom. In the more general case (12.661) it is 
practically certain that imposing En ^ constrains the system too much. A somewhat 
similar problem has been discovered in the original Hofava-Lifshitz theory jl3j . 

3 Conclusion 

As a summary, there are serious problems with the physical consistency of the non- 
projectable version of our theory, which will likely be impossible to resolve with addi- 
tional constraints. For the general potential f l2.62p the conclusion is similar compared 
to the usual Hofava-Lifshitz gravity: since iV = is required by the preservation of the 
Hamiltonian constraints under time evolution, the theory is physically inconsistent. The 
difference is that in our E{R) gravity version, the undesirable condition = cannot be 
avoided consistently even for the low-energy effective action. In the light of Ref . [26] , it 
is very unlikely that there would exist a form of the action with scaling dimension z = 3 
that could avoid this problem in the usual Hofava-Lifshitz gravity. The same is expected 
to apply in our theory. 

Thus, we can conclude that only the version of the modified F{R) Hofava-Lifshitz 
gravity with a projectable lapse, A^ = N(t), is a practicable theory of modified gravity. 
The non-projectable version of the theory is troubled by a similar physical inconsistency 
as the usual non-projectable Hofava-Lifshitz gravity. 

It would be very interesting to understand whether some further generalization of the 
Hofava-Lifshitz gravity could be defined consistently without the projectability condition 
on the lapse — perhaps even along the lines of the general modified first-order Hofava- 
Lifshitz gravity of Ref. [17] . It would require that the coefficients E^ in f l2.64p -f l2.66p take 
a very special form, preferably vanish altogether. Another intriguing prospect would be 
a generalization of the "healthy extension" of usual Hofava-Lifshitz gravity [10] to our 
modified F[R) Hofava-Lifshitz gravity. 
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